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We study the Hartree ground state of a dipolar condensate of atoms or molecules in an three- 
dimensional anisotropic geometry and at T = 0. We determine the stability of the condensate as a 
function of the aspect ratios of the trap frequencies and of the dipolar strength. We find numerically 
a rich phase space structure characterized by various structures of the ground-state density profile. 



I. INTRODUCTION 

The realization of Bose-Einstein condensation (BEC) 
in low-density atomic vapors [l|, HJ has led to an ex- 
plosion of experimental and theoretical research on the 
physics of quantum-degenerate atomic and molecular sys- 
tems. While much of the work so far has concentrated on 
systems characterized by s-wave two-body interactions, 
the recent demonstration of a condensate of chromium 
atoms [1] opens up the study of gases that interact via 
long-range, anisotropic magnetic dipole interactions. In 
a parallel development, it can be expected that quan- 
tum degenerate samples of heteronuclear polar molecules 
will soon be available through the use of Fcshbach reso- 
nances 0, [1[ , photoassociation [f| 0] , or a combination of 
the two. When in their vibrational ground state, these 
molecules interact primarily via the electric dipole poten- 
tial, and they are expected to provide a fascinating new 
type dipole-dominated condensates in the near future. 

As a result of the anisotropy and long-range nature of 
the dipole-dipole interaction, a number of novel phenom- 
ena have been predicted to occur in low-density quantum- 
degenerate dipolar atomic and molecular systems, both 
in conventional traps and in optical lattices. An early 
study of the ground state of polar condensates was pre- 
sented in Ref . [§[ , which determined its stability diagram 
as a function of the number of atoms and s-wave scatter- 
ing length. It identified a stable structured ground state 
for a specific range of parameters. At about the same 
time, the effect of trap geometry on the stability of the 
condensate was considered in Ref. @ for a system domi- 
nated by the dipole interaction. This was followed by the 
prediction [irj| of the existence of a number of quantum 
phases for dipolar bosons in optical lattices. Recent work 
[111 . [T3 | considers the structural phases of vortex lattices 
in rotating dipolar Bose gases. 

A novel feature of dipolar condensates, as compared to 
their scalar cousins, is the appearance of a roton mini- 
mum in their Bogoliubov spectrum. This feature was dis- 
cussed in the context of atomic condensates in Ref. [HI], 
which considered the impact of the roton-maxon feature 
in the excitation spectrum and the stability of pancake- 
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shaped dipolar condensates. For this particular geometry 
it was found that the excitation spectrum can touch the 
zero-energy axis for a non-zero wave vector which 
points to the instability of homogeneous condensates and 
the onset of density modulations [To| . A roton minimum 
was also found |16| for the case of laser-induced dipolar 
interactions in self-bound BECs with cylindrical symme- 
try. Quasi-2D dipolar bosons with a density-modulated 
order parameter were determined to be unstable within 
the mean-field theory [TtJ , and cigar-shaped quasi-one di- 
mensional condensates were likewise found [18| to be dy- 
namically unstable for dipolcs polarized along the axis of 
the cylindrical trap. The stability of dipolar condensates 
in pancake traps was also recently discussed in Ref. [2l| , 
which found the appearance of biconcave condensates for 
certain values of the trap aspect ratio and strength of the 
dipole interaction. From the Bogoliubov excitation spec- 
trum it was possible to attribute the instability of the 
condensate under a broad range of conditions to its az- 
imuthal component. 

Further building on these studies, the present note 
reports the results of a detailed numerical analysis of 
the stability and structure of the Hartree ground state 
of dipolar condensates confined in anisotropic harmonic 
trap. We proceed by introducing the trap frequencies 
u> x ,LUy and u) z , respectively, in the x, y and z directions, 
and the corresponding trap aspect ratios X y = ui y /ui x 
and A z = lu z /lu x . Thus X z = 1 corresponds to a pancake 
trap, whereas X z — corresponds to a cylindrical trap 
with free motion in z direction. We further assume for 
concreteness that an external field polarizes the dipolcs 
along the y axis. The stability of the condensate is then 
determined as a function of the trap aspect ratios and 
of an effective dipolar interaction strength that is pro- 
portional to the number of atoms or molecules in the 
condensate. Various ground state structures of the con- 
densate are identified in the stable region of parameter 
space. 

The remainder of this paper is organized as follows: 
Section II introduces our model and comments on im- 
portant aspects of our numerical approach. Section III 
summarizes our results, identifying up to five different 
types of possible ground states, depending on the tight- 
ness of the trap and the particle number. Finally, Section 
IV is a summary and conclusion. 
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II. FORMAL DEVELOPMENT 

The dipole-dipole interaction between two particles 
separated by a distance r is 



V dd {r) = .9dd 



1 - 3y 2 /r 2 



(1) 



where <?dd is the dipole-dipole interaction strength and y 
is the polarization direction. For atoms with a permanent 
magnetic dipole moment we have <7dd = Mo/-* 2 n /47r while 
for dipolar molecules gdd = A t e/4 7re o J Mm and being 
the magnetic moment of the atoms and the electric dipole 
moment of the molecules, respectively. 

Within the mean-field approximation, the condensate 
order parameter <fi(r) satisfies the Gross-Pitaevskii (GP) 
equation 



E<j>(r) = [H +g\^(r)\ 

+ N f V dd (r-r')\^(r')fd 3 r' 



where 



2m 



0(0, (2) 



V 2 + V(, 2 + Ay + A^) (3) 



is the sum of the kinetic energy and the trapping poten- 
tial and TV denotes the number of particles in the conden- 
sate. The second term on the right-hand side of Eq. @ 
is the contact interaction, g = 47rft 2 a/m being propor- 
tional to the s-wave scattering strength a, and the third 
term describes the effects of the nonlocal dipole-dipole 
interaction. For dipole interaction dominated systems, 
g is small compared to V dd . This is the case that we 
consider here, and in the following we neglect the s-wave 
scattering term altogether. 

For convenience we introduce the dimensionlcss param- 
eter 



D = Ng dd m/{£ X H 2 ) 



(4) 



that measures the effective strength of the dipole-dipole 
interaction, where the oscillator length £ x = yj h/{muj x ). 
The condensate ground state is then determined numer- 
ically by solving the Gross-Pitaevskii equation © for 
imaginary times. The term involving the dipole interac- 
tion energy is calculated using the convolution theorem, 



V dd (r-r')\cj)(r')\ 2 d 3 r' = 



where T and T~ x stand for Fourier transform and in- 
verse Fourier transform, respectively. The dipole-dipole 
interaction is calculated analytically in momentum space 
as B, 



r[V dd (r)\ 



47T 

T 



V2 



■*2 



k 2 



(5) 



k x , k y , k z are the momentum components in x, y, z direc- 
tion. 

The initial order parameter was chosen randomly, and 
the stability diagram was generated for each pair of 
parameters (Aj,,A 2 ) by increasing the effective dipolar 
strength D until a critical value D cr above which the 
condensate collapses. Because of the random initial con- 
dition this value varies slightly from run to run. The plot- 
ted results show the average over 100 realizations of the 
initial wave function, the error bars indicating the maxi- 
mum deviation from of D cl - from its mean D cr . This ap- 
proach typically resulted in numerical uncertainties sim- 
ilar to those of Ref . . 



III. RESULTS 

A good starting point for the discussion of our results 
is the observation that in the case of a cylindrical trap, 
A 2 = 0, wc found no stable structured condensate ground 
state. (By structured profiles, we mean profiles that are 
not simple gaussians.) In particular, solutions exhibit- 
ing density modulations along the z-axis were found to 
be unstable. Moving then to the case of a pancake trap 
by keeping X y fixed but increasing A 2 from to 1, we 
found for X y > 4 the appearance of a small parameter re- 
gion where the stable ground state is characterized by a 
structured density profile, the domain of stability of this 
structured solution increasing with X y . A {A 2 — D} phase 
space stability diagram typical of this regime is shown in 
Fig. [1] for X y = 5. In this figure, region / is characterized 
by the existence of a usual condensate with its familiar, 
structureless gaussian-like density profile. As A 2 is in- 
creased, the condensate becomes unstable for decreasing 
values of the effective dipole interaction strength D, or al- 
ternatively of the particle number TV. For .525 < A 2 < .7, 
though, the ground state changes from a gaussian-shaped 
to a double-peaked density profile (domain II in the fig- 
ure), before the system becomes unstable. 

Figures 0] and O show surface plots and corresponding 
3-D renditions of density profiles typical of the various 
situations encountered in our study. Figures [4^, and [5^ 
are illustrative of the present case. The appearance of 
two density peaks away from the center of the trap results 
from the interplay between the repulsive nature of the 
dipoles in a plane transverse to its polarization direction, 
the (a; — z) plane, and the confining potential. 

Increasing the tightness of the trap along the polariza- 
tion direction y, i.c, increasing X y , results in the emer- 
gence of additional types of structured ground states. 
One such case is illustrated in Fig. [2l which is for A 2 = 
5.5. For small values of A 2 , i.e., a weak trapping po- 
tential along the z-direction, we observe the appearance 
of a domain (region III in the figure) characterized by 
a double-peaked ground state with the maximum den- 
sity along the z direction and a gaussian-likc density in 
x direction. This type of double-peaked structure along 
the weak trapping axis was first predicted in Ref. [20| . 
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FIG. 1: (\ Z ,D) stability diagram of a dipolar condensate in 
an anisotropic trap for X y = 5. The condensate is unstable 
in the region above the solid line. The dashed line is the 
boundary between a "structureless gaussian" and a double- 
peaked ground-state density profile. The error bars give an 
indication of the accuracy of the numerical simulations. 



Typical density profiles in this region resemble those in 
Figs. 4a and 5a, but with a rotation by 90 degrees in the 
(x, z)-planc. The regions II and III are separated by a 
small additional domain IV characterized by a ground- 
state distribution with a quadruple-peaked structure as 
illustrated in Figs. 4b and 5b, as might be expected. In 
general, these characteristics of the ground state density 
profile persist until X y ~ 6.5. 

Figure [3] shows a stability diagram typical of higher 
values of the aspect ratio X y , in this case X y = 7, for 
0.4 < X z < 1. As D is increased, the ground-state den- 
sity of the condensate first undergoes a transition from 
a gaussian-like to a double-peaked profile of the type il- 
lustrated in Figs. 4b and 5b (region III). As D is further 
increased, this domain is followed for X z close enough to 
unity by a second transition to a domain (region V) with 
the appearance of a density minimum near trap center. 
Initially, this minimum is surrounded by a region with a 
radial density modulation, see Fig. 4c, but for larger val- 
ues of X z this modulation is reduced, see Figs. 4d and 5d. 
In that region, the density profile resembles the solution 
previously reported in Rcf. plj for a similar parameter 
range. 

In the case of atoms a typical magnetic moment of 6[1b, 
and we find that the range of critical dipolc strengths 
D corresponding to structured ground states can be 
achieved for 10 4 — 10 5 atoms for trapping frequencies 
cu x w 1kHz. For molecules with a typical electric dipole 
moment of IDebye the corresponding number is 10 3 — 10 4 
molecules. While these are relatively high particle num- 
bers, especially for the atomic case, they do not seem out 
of reach of experimental realization. 
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FIG. 2: Stability diagram of a dipolar condensate in an 
anisotropic trap with \ y — 5.5, as a function of the dipolar 
strength D and aspect ratio X z . The region above the solid 
line corresponds to unstable solutions of Eq.Q. The dashed 
line is the boundary between the structured and the standard 
ground-state density profiles. The dot-dashed line marks the 
boundary between the domains with quadruple-peaked and 
double-peaked ground-state density profiles. 



IV. CONCLUSION 



In conclusion, we have performed a detailed numerical 
study of the ground state structure and stability of ul- 
tracold dipolar bosons in an anisotropic trap for dipolcs 
polarized along the ^-direction. The trap aspect ratios 
along y and z direction, X v and X z were used as control 
parameters, and the mean-field stability diagram has es- 
tablished as function of these parameters and a dimen- 
sionlcss interaction strength D. For small X y the system 
was found to exhibit a standard density profile, but for 
larger values, and depending on X z , various structured 
ground state were found to appear before reaching the 
unstable regime where the condensate collapses. These 
include a four-peak structured solution in the x — z plane, 
a ring-like ground state with a modulated radial density 
profile. For X y ~ 7 and X z = 1, we found a biconcave 
condensate profile, as already reported in [2l| . 

For strong confining potentials along the dipolc polar- 
ization direction, i.e for large X y , increasing X z can be 
viewed as resulting in a change from a quasi-one dimen- 
sional to a quasi-two-dimensional geometry. As such we 
can think of the various ground-state structures as a re- 
sult of dimensional crossover in a trapping geometry. To 
gain a deeper understanding of these structures as we 
approach the instability region, future work study the 
Bogoliubov spectrum of the trapped system. 
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FIG. 3: Stability diagram of a dipolar condensate in an 
anisotropic trap with X y — 7, as a function of the dipolar 
strength D and the aspect ratio \ z . The region above the 
solid line is characterized by unstable ground-state solutions 
of the mean-field equation The dashed line marks the 

boundary between a structured and a gaussian-like ground- 
state density profile. The dotted line is the boundary between 
the region with ring-like and two-peaked condensate in (x — z) 
plane. 
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FIG. 4: Two-dimensional surface plots of the structured ground-state density profiles typical of various stable domains: (a) 



Double-peaked density profile in the 



z)-plane for region II, for the parameters X y — 5, A z = .6 and D — 23. The points 



of maximum density are away from the trap center and along the a;-direction. In region (III) the shape of the condensate is 
similar, but with maximum density along z direction, (b) Typical quadruple-peaked density profile characteristic of region IV. 
Here X y — 5.5, A z = .575 and D = 42 (c) Stable ground state solution in region V with \ z < 1. The density is higher and 
modulated on a radius away from the trap center in the x — z plane. In this example Aj, = 7, Az = .85 and D = 40. (d) 
Density profile typical of the domain V. The condensate is biconcave with maximum density along a constant radius from the 
trap center. In this example, X y — 7 and D = 32. 




FIG. 5: Three-dimensional ground-state density profile for the same parameters an in Fig. [4] 



